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Abstract 

The differential cross sections for the process 77 — > 7r°7r° have been measured in the kinematic 
range 0.6 GeV < W < 4.1 GeV, | cos(9*| < 0.8 in energy and pion scattering angle, respectively, in 
the 77 center-of-mass system. The results are based on a 223 fb _1 data sample collected with the 
Belle detector at the KEKB e + e~ collider. The differential cross sections are fitted in the energy 
region 1.7 GeV < W < 2.5 GeV to confirm the two-photon production of two pions in the G 
wave. In the higher energy region, we observe production of the XcO charmonium state and obtain 
the product of its two-photon decay width and branching fraction to 7r°7r°. We also compare the 
observed angular dependence and ratios of cross sections for neutral-pion and charged-pion pair 
production to QCD models. The energy and angular dependence above 3.1 GeV are compatible 
with those measured in the ir + Tr~ channel, and in addition we find that the cross section ratio, 
o-(7T 7r°)/V(7r + 7r-), is 0.32 ± 0.03 ± 0.05 on average in the 3.1-4.1 GeV region. 

PACS numbers: 13.20.Gd. 13.60.Le, 13.66.Bc, 14.40. Cs,14.40.Gx 
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I. INTRODUCTION 



Measurements of exclusive hadronic final states in two-photon collisions provide valuable 
information concerning the physics of light and heavy-quark resonances, perturbative and 
nonperturbative QCD, and hadron-production mechanisms. So far, we have measured the 
production cross sections for charged-pion pairs [l|, 0] , charged- and neutral- kaon pairs @, 
y, Ej], an d proton-antiproton pairs jEJ. We have also analyzed D-meson-pair production 
and observe a new charmonium state 0. Recently, we have presented a measurement 
of neutral-pion pair production based on a data sample corresponding to an integrated 
luminosity of 95 fb _1 [7|. We have carried out an analysis in the energy range W < 1.6 GeV 
to extract information on light quark resonances from the energy and angular dependence 
of the differential cross sections, by fitting to the resonance parameters of the /o(980), 
/2(1270) and additional hypothetical resonances. The statistics of these measurements is 2 
to 3 orders of magnitude higher than in the pre-S-factory measurements opening a new 
era in studies of two-photon physics. 

Here we present measurements of the differential cross sections, da / d\ cos 9* \, for the 
process 77 — > 7r 7r° in a wide two-photon center-of-mass (cm.) energy (W) range from 0.6 
to 4.1 GeV, and in the cm. angular range, | cos^*] < 0.8. We use a 223 fb _1 data sample, 
which is more than twice as large as that in our previous analysis [t|. We focus on the range 
W > 1.4 GeV, where the previous data was statistically limited. 

In the intermediate energy range (1.0 GeV < W < 2.4 GeV), production of two pions is 
dominated by intermediate resonances. For ordinary qq mesons in isospin conserving decays 
to 7r7r, the only allowed I G J PC states produced by two photons are + (even) ++ , that is, 
fj =even mesons. Several mesons with these quantum numbers are suggested by results of 
hadron-beam or charmonium decay experiments in the 1.5 - 2.2 GeV region. However, none 
of them have been firmly established in two-photon processes, which are sensitive to the 
internal quark structure of the meson. In addition, the 7r°7r° channel has two advantages in 
the study of resonances: a smaller contribution from the continuum is expected in it than 
in the 7r + 7r~ channel; and the angular coverage is larger (| cos#*| < 0.8 instead of 0.6). 

At higher e nerg ies (W > 2.4 GeV), we can invoke a quark model. In leading-order 
calculations 0, [To[ 11], which take into account the spin correlation between quarks, the 



7r°7r° cross section is predicted to be much smaller than that of tt + tt~, suggesting a ratio 
of 7r°7r° to 7r + 7r~ cross sections around 0.04-0.07. However, higher-order or nonperturbative 
QCD effects can modify this prediction. For example, the handbag model, which considers 
soft hadron exchange, predicts the same amplitude for the two processes, and thus the 
expected ratio is 0.5 [l2j. Analyses of energy and angular distributions of the cross sections 
are essential for determining the properties of the observed resonances and for testing the 
validity of QCD models. 

The organization of this article is as follows. In Sec. [TTJ a brief description of the Belle 
detector is given. Section II III explains the procedure used to obtain differential cross sec- 
tions. Section IIVI is devoted to results on the two-photon production of two pions in the G 
wave obtained by fitting differential cross sections in the range 1.7 GeV < W < 2.5 GeV. 
Section [V] describes analyses at higher energy. The topics included there are the angular 
dependence as a function of W, the observation of the Xco and Xc2 charmonia states and the 
ratio of cross sections for 7r°7r° to ti + tx~ production. Finally, Sec. I VI I summarizes the results 
and presents the conclusion of this paper. 
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TABLE I: Data sample: luminosities and energies 



e + e cm. energy 


Luminosity 


Comment 


fGeV) 






10.58 


179 


T(4S) runs 


10.52 


19 


continuum runs 


10.36 


2.9 


T(3S) runs 


10.30 


0.3 


continuum runs 


10.86 


21.7 


T(5S) runs 


total 


223 





II. EXPERIMENTAL APPARATUS 



We use a 223 fb 1 data sample from the Belle experiment [13J at the KEKB asymmetric 



energy e ' e collider [14|. The data were recorded at several e + e~ cm. energies summarized 
in Table [H The difference of the luminosity functions (two-photon flux per e + e~-beam 
luminosity) in the measured W regions due to the difference of the beam energies is small 
(maximum ± 4%). We combine the results from the different beam energies. The effect on 
the cross section is less than 0.5%. 



The analysis is carried out in the "zero-tag" mode, where neither the recoil electron nor 
positron are detected. We restrict the virtuality of the incident photons to be small by 
imposing strict transverse-momentum balance with respect to the beam axis for the final- 
state hadronic system. 



A comprehensive description of the Belle detector is given elsewhere [13| . We mention here 



only those detector components that are essential for the present measurement. Charged 
tracks are reconstructed from hit information in the silicon vertex detector and the central 
drift chamber located in a uniform 1.5 T solenoidal magnetic field. The detector solenoid is 
oriented along the z axis, which points in the direction opposite to that of the positron beam. 
Photon detection and energy measurements are performed with a CsI(Tl) electromagnetic 
calorimeter (ECL). 

For this all-neutral final state, we require that there be no reconstructed tracks coming 
from the vicinity of the nominal collision point. Therefore, the central drift chamber is used 
for vetoing events with charged track(s). The photons from decays of two neutral pions are 
detected and their momentum vectors are measured by the ECL. The ECL is also used to 
trigger signal events. 



III. DERIVING DIFFERENTIAL CROSS SECTIONS 

The event triggers, data processing, and event selection are the same as those described 
in Ref. Q. We derive the cm. energy W of the two-photon collision from the invariant mass 
of the two-neutral-pion system. We calculate the cosine, | cos 9* | of the 7r° scattering angle 
in the 77 cm. frame for each event, using the e + e~ collision axis in the e + e~ cm. frame 
as the reference axis for the polar angle. The possible bias due to the unknown 77 collision 
axis is negligible. 
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A. Data reduction 



We find that the signal candidates in the low energy region (W < 1.2 GeV) are con- 
siderably contaminated by background. In order to separate the signal and background 
components, we study the p t -balance distribution, i.e., the event distribution in | J2Pt I- We 
estimate the p t - unbalanced background component for W < 1.2 GeV in the same manner 
as in the previous analysis 0] and subtract the yield in the signal region. However, above 
1.2 GeV, we cannot quantitatively determine the background contamination because of the 
small background rate and low statistics of the sample, as well as the uncertainty in the 
functional form for the signal shape. 

Using the ratio of yields between the pt-balanced and unbalanced regions, we can estimate 
the backgrounds. In Fig. [TJ we plot the W dependence of R defined as: 

F(0.15 GeV/c < | £p?| < 0.20 GeV/c) 

y(|EP*| <0.05 GeV/c) ' U 

where Y is the yield in the indicated | Pt I region. We integrate over all angles in this 
figure. The main part of the W dependence of R comes from the energy dependence of the 
momentum resolution. The expected ratio from the pure signal component is shown by the 
solid line. The signal events for e + e~ —>■ e + e~ir ir are generated using the TREPS code (IBj . 
All Monte Carlo (MC) events are put through the trigger and detector simulators and the 
event selection program. The MC events are corrected for MC/data difference in the p t 
resolution discussed in the next section. The excess of R over the line (AJf2) is expected to 
correspond to the contribution from the pj-unbalanced background. The excess is relatively 
small above 1.0 GeV, although some fine structure is visible there. In the range 1.2-3.3 GeV, 
AR ranges between 0.00 and 0.08, and above 3.3 GeV it is in the range from 0.08 to 0.2. 
From the R values, we estimate that the background contamination in the signal region is 
~ -R/4, which is smaller than 3% for 1.5 - 3.3 GeV and around 3% for 3.6 - 4.1 GeV. We 
subtract 3% for 3.6 - 4.1 GeV, and assign a 3% systematic error from this source for the full 
1.5 - 4.1 GeV range. 

We estimate the invariant-mass resolution from studies of signal-MC and experimental 
distribution. The true W distribution in the range 0.9 GeV < W < 2.4 GeV is obtained 
by unfolding the differential cross sections as described in Ref. 0]. For lower energies, 
W < 0.9 GeV, the effect of the migration is expected to be small because the invariant-mass 
resolution is much smaller than the bin width. For higher energies, W > 2.4 GeV, where the 
statistics is relatively low and unfolding would enlarge the errors, we adopt a rather wide 
bin width (100 MeV) without unfolding. A total of 2.90 x 10 6 events are selected in the 
region of 0.6 GeV < W < 4.1 GeV and | cos#*| < 0.8. 



B. Calculation of differential cross sections 

We determine the efficiency for the si gna l using a full MC simulation. The MC signal 



events generated using the TREPS code [15j are isotropically distributed in | cos6*\ at 58 



fixed W points between 0.5 and 4.5 GeV. The angular distribution at the generator level 
does not play a role in the efficiency determination, because we calculate the efficiencies 
separately in each | cos#*| bin with a 0.05 width. Samples of 4 x 10 5 events are generated 
at each W point. Two sets of different background conditions, which were extracted from 
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FIG. 1: The yield ratio R in the pt-unbalanced bin to the pt-balanced (signal) bin (see text 
for the exact definition) for the experimental data. The solid line shows the signal component 
obtained from the signal MC and corrected taking into account the poorer momentum resolution 
in experimental data. 

the beam collision data are embedded in the signal-MC data in the detector simulation. To 
minimize statistical fluctuations in the MC calculation, we fit the numerical results of the 
trigger efficiency to a two-dimensional empirical function in (W, | cos#*|). 

The efficiency calculated from the signal-MC events is corrected for a systematic difference 
of the peak widths in the p r balance distributions found between the experimental data and 
the MC events, which is attributed to a difference in the momentum resolution for 7r°'s. The 
correction factor is typically 0.95. 

The differential cross section for each (W, | cos#*|) point is given by: 

da AY -AB 

d\cos6*\ ~ AW A\ cos 6*\j CdtL^(W)r! ' ^ ' 

where AY and AB are the signal yield and the estimated p t -unbalanced background in the 
bin, AW and A| cos^*! are the bin widths, / Cdt and L^(W) are the integrated luminosity 
and two-photon luminosity function calculated by TREPS [15], respectively, and i] is the 
efficiency including the correction described above. The bin sizes for W and A| cos#*| are 
summarized in Table HT1 



Figure [2] shows the W dependence of the cross section integrated over | cos6**| < 0.8. We 
have removed the bins in the range 3.3 GeV < W < 3.6 GeV, because we cannot separate 
the Xco an d Xc2 components and the continuum in a model-independent way due to the 
finite mass resolution and insufficient statistics of the measurement. The cross section in 
this region is discussed in detail in Sec. [Vj 
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TABLE II: Bin sizes 



W range 


AW 


A| cos 0*| 


(GeV) 


(GeV) 




0.6 - 1.8 


0.02 


0.05 


1.8 - 2.4 


0.04 


0.05 


2.4 - 4.1 


0.10 


0.05 



1000 



100 



10 



o 

CD 
O 



0.1 



0.01 



J\ 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 

W (GeV) 



FIG. 2: The integrated cross section in the angular regions | cos 6*\ < 0.8. Data points in bins near 
3.5 GeV are not shown because of uncertainty from the Xcj subtraction. 



We show the angular dependence of the differential cross sections at several W points in 
Fig. [3l Note that the cross sections in neighboring bins after the unfolding are no longer 
independent of each other in either central values or size of errors. 

The systematic uncertainties for the cross sections arise from various sources; they are 
listed in Table III II together with the estimated values. Uncertainties from the unfolding 



procedure, using the single value decomposition approach in Ref. [16j, are estimated by 
varying the effective-rank parameter of the decomposition within reasonable bounds. 

The total systematic error is obtained by adding the uncertainties in quadrature and is 
about 10% in the intermediate W region (1.04 GeV < W < 3.0 GeV). It becomes much 
larger at lower W . At higher W, the systematic error is rather stable, typically about 11%. 



IV. STUDY OF G-WAVE ACTIVITY 

Previously, we have obtained a reasonable fit to a simple model of resonances and smooth 
backgrounds in the energy region 0.8 GeV < W < 1.6 GeV from the differential cross 
sections of 77 — > 7r°7r° with a 95 fb _1 data sample Q. The clear / (980) peak and the 
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W = 0.97 GeV 
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0.2 0.4 0.6 O.f 





W= 3.15 GeV 



If 



i'ff 



ill 



0.2 0.4 0.6 
| cos B*| 



0.4 0.6 

cos e*| 



FIG. 3: The differential cross sections for W points indicated. The bin sizes are summarized in 
Table IE 

TABLE III: Systematic errors for the differential cross sections. Ranges of errors are shown when 
they depend on W. 



Source 



Error (%) 



Trigger efficiency 4-30 
7T° reconstruction efficiency 6 

Pt-balance cut 1.5-5 

Background subtraction 0-40 

Luminosity function 4-5 
Overlapping hits from beam background 2-4 
Other efficiency errors 4 

Unfolding procedure 0-4 



Overall 



typ. 10 - 11 



large contribution from the /2(1270) can be fitted with parameters determined from 7r + 7r~ 
data [l|. 

In this section, we concentrate on the G wave, in particular, the ^4(2050) resonance, 
whose existence is well established, but whose production in two-photon collisions has never 
been positively identified. We fit the energy region 1.7 GeV < W < 2.5 GeV using a 
high-statistics sample of 223 fb" 1 that contains 2.3 times more events than in the previous 
experiment Q; the number of events in this region is 155k. 



A. Parametrization of Partial Wave Amplitudes 

In the energy region W < 3 GeV, J > 4 partial waves (the next is J = 6) may be 
neglected so that only S, D and G waves are to be considered. The differential cross section 
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FIG. 4: Spectrum of 5 (top section of left-hand panel), -Do (bottom section of left-hand panel) 

* 2 

and L>2 (right panel) for 1.7 GeV < < 2.5 GeV and results of parametrization (see text). The 
error bars shown are diagonal statistical errors. 



can be expressed as 
da 



dQ 



(77 -> ttV) 



SY» + D Y» 



D 2 Y 2 2 



G 2 y} 



(3) 



where Dq and Go {D 2 and G2) denote the helicity (2) components of the D and G waves, 
respectively, and Y™ are the spherical harmonics. Since the |lj"|s are not independent, 
partial waves cannot be separated using measurements of differential cross sections alone. 
To overcome this problem, we write Eq. ([3]) as 



da 



-(77 -> vrV) 



Aird\ cos 9* 



S 2 \Y °\ 2 



f)2 |y0|2 



Di\Y 2 r +G 2 \Yl 



1 2 



A2 \\r2 1 2 

u- 2 |r 4 1 



(4) 



The amplitudes S 12 , etc. correspond to the cases where interference terms are neglected; 
they can be expressed in terms of S, D , D 2 , Go, and G 2 0). Since squares of spherical 
harmonics are independent of one another, we can fit differential cross sections at each W 
to obtain S* 2 , Z)q, D\, G\, and G\. For | cos6*\ < 0.7, the |F 4 °| 2 and \Y 2 \ 2 terms are nearly 
equal, so we fit Gq + G\ and G\ — G\ instead. The resulting spectra are shown in Figs. 0] 
and [3 

We parametrize the partial wave amplitudes in terms of resonances and smooth "back- 
grounds". Once the functional forms of the amplitudes are fixed, we can use Eq. ()3]) to 
fit the differential cross sections. From Fig. [51 it appears that the G wave contributions 
are nonzero for W ~ 1.8 GeV and are dominated by the G2 wave. Here we assume (and 
check the necessity of) including the ^4(2050) in the G2 wave. Since the G2 wave interferes 
with the D 2 wave, we include the resonance /2(1950), which is known to couple to two pho- 
tons [!, [l7[ • There are several other resonances that might couple to 77 and irir in this mass 
region, which are listed in Ref. 17|. Here we assume that the /2(1950) is just an empirical 
parametrization representing these other resonances; we denote it here as the "/2(1950)". 
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FIG. 5: Spectrum of Go + G2 (top section) and Go — G2 (bottom section) for 1.7 GeV < W < 
2.5 GeV and fitted curves (see text). The error bars shown are diagonal statistical errors. 

We parametrize the partial waves as follows: 



5 


= B s , 


D 


— B D0 , 


D 2 


= A "fa (1950) 


Gq 


= 0, 


G 2 


= ^4/ 4 (2050)e 



»e J<p2 +B 



D2, 



+- B G2 , (5) 

where ^4"j 2 ( 1950 )" an d ^4/ 4 (2050) are the amplitudes of the corresponding resonances; Bs, -Bdo, 
B02 and Bq2 are nonresonant (background) amplitudes for S, Do, D2 and G2 waves; and 
<p2 and 04 are the phases of resonances relative to background amplitudes. We assume that 
Go = and that G2 consists only of the ^(2050) and a smooth background. 

The relativistic Breit-Wigner resonance amplitude Ar{W) for a spin-J resonance R of 
mass rriR is given by 



A J R (W) 



'8vr(2J+ l)m R 



r tot r 77 i3(7r°7r°) 



w 



m. 



— W 2 — imjiT 



(6) 



tot 



The resonance parameters given in Ref. [17j for the ^2(1950) and ^(2050) are summarized in 
Table HVl We assume an energy-independent width for the "/ 2 (1950)" and / 4 (2050) because 
most of their individual decay fractions are unknown. 



The background amplitudes are parametrized as follows. 

B s = a sr (W - Wo) 2 + b sr (W - W ) + c sr + % (a si (W - W ) 2 + b si (W - W ) + c si ) , 
B D0 = a (W - W ) 2 + b (W - W ) + c , 

B m = a 2r (W - Wo) 2 + b 2r (W - W ) + c 2r + 1 (a 2l {W - W ) 2 + b 2l (W - W ) + c 2i ) , 
B G2 = a g (W - Wo) 2 + b g (W - Wo) + c g (7) 
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TABLE IV: Parameters of the / 2 (1950) and .A (2050) 0]. 



Parameter / 2 (1950) / 4 (2050) Unit 



Mass 


IMA ± 12 


2018 ± 11 


MeV/c 


Width 


472 ± 18 


237 ± 18 


MeV 


B(tttt) 


seen 


17.0 ±1.5 


% 


B{KK) 


seen 


0-68±8;?| 


% 




seen 


0.21 ±0.08 


% 




seen 


unknown 





where Wq = 1.7 GeV. The background amplitudes D and G 2 are taken to be real by 
definition. The other background amplitudes are assumed to be quadratic in W for both 
their real and imaginary parts. We fix Bq2 — at W — 1.7 GeV (c g = 0) to reduce the 
number of parameters; leaving c g free does not improve the fits. 



B. Fit results 

We minimize y 2 defined as 

gr^yCWi, I COSnOdata " ar^W, I COS0*|,-) pred . 



x 2 = E 



dir 



-AW h | COS0*|j)data 



i,j \ d\cos9 

where the summation is over (Wi, | cos#*|j) bins, da/d\ cos9*\(Wi, | cos#*|j)data (pred.) is the 
cross section data (prediction using Eq. (l5"]) ) at a bin (Wi, | cos#*|j), and the denominator is 
the estimated statistical error. 

When the mass and width of the ^4(2050) are fixed to the values given in the PDG ta- 



bles 17J] as summarized in Table HV] then the fit is very poor yielding \ 2 (ndf) = 594.4 (313) 
(see Table IVl) . This is to be compared with 323.2 (311) obtained when the mass and width 
are floated. In this paper we quote the results of the fits with the mass and width of the 
/ 4 (2050) as free parameters. 

Here the unfolded differential cross sections are fitted. Fits are performed 1000 times 
for each study with randomly-generated initial values for the parameters A unique solution 
with good quality (x 2 / n df = 1.04) is repeatedly found ("nominal fit"). The fit results are 
shown in Fig. [6] for the differential cross sections, in Fig. [7] for the total cross section and in 
Figs. H] and [5] for the S 2 , etc. Since the two-photon coupling of the ^(2050) has not been 
measured before, a fit without this resonance is also given in Table |V] The fit quality is 
unacceptable, strongly indicating that the ^(2050) has a nonzero two-photon coupling. A 
fit without the "/2(1950)" is also made giving a much worse fit and is included in Table IVl 

We have performed additional fits to investigate whether we can conclude that the 
^(2050) is mainly produced in the helicity-2 state. Note that the angular dependence 
of F 4 ° and F 4 2 is very similar for | cos6*\ < 0.7 and hence it is expected to be rather difficult 
to distinguish Go and G 2 waves. A fit where the role of Go and G 2 is interchanged (i.e. by 
setting G2 = and by including the ^(2050) and background in Go) yields \ 2 — 448.2, 
which can be compared to 323.2 for the nominal fit. However, more reasonable fits are 
obtained when two more parameters are introduced in the G 2 background (a second-order 
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FIG. 6: Measurements of da/d\ cos 9*\ (nb) (data points) and results of the fit (thick solid line) 
for the VF-bins indicated. The dotted, dashed, dot-dashed and thin lines indicate \S\ 2 , A7r\DoY 2 \ 2 , 
4ir\D 2 Y 2 2 \ 2 , 4tt|G 2 14 2 | 2 , respectively. 



polynomial real amplitudes or linear amplitudes for both real and imaginary parts, which 
are set to zero at W = 1.7 GeV). Here the fits where the role of G 2 and Go have been 
interchanged are also compared. In each case, the ^(2050) in a G 2 wave is favored over 
that in a Go wave with a x 2 differences of about 6 and 26, respectively. Thus, the helicity-2 
production of the ^(2050) is favored but not conclusively. 

According to Fig. [71 the \D 2 \ 2 term has an enhancement around W = 2.35 GeV, which 
might be identified as the f 2 (2300). To study this possibility, a fit is made including the 
f 2 (2300). The conclusion is, however, that we have no sensitivity to it; the fit does not 
improve significantly by its inclusion. We believe that the enhancement arises from the 
"/2(1950)" and its interference with the G2 wave and underlying continuum, i.e., a fit without 
the "/2(1950)" gives a smooth D 2 amplitude (with much worse \ 2 (Table IV]) ). 



C. Study of systematic errors 

Various sources of systematic errors on the parameters are considered such as dependence 
on the fitted region, normalization errors of the differential cross sections, assumptions on 
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FIG. 7: Total cross section (| cos#*| < 0.8) (nb)) and the results of the nominal fit (curves). 



TABLE V: Fitted parameters 



Parameter 


Nominal 


Fixed / 4 (2050) No / 4 (2050) No 7 2 (1950)" 


Unit 


Mass(/ 4 (2050)) 


1885lJ| 


2025 (fixed) 




2052 ± 6 


MeV/c 2 


r t ot(/ 4 (2050)) 


453 ± 20 


225 (fixed) 




2571? 


MeV 


r 77 i3(7r°vr ) 


7 7+1.2 
'■'-1.1 


11.8 ±0.2 


(fixed) 


14 2 +0 - 9 


eV 


Mass(7 2 (1950)") 


20381}? 


20261 2 


21141 J| 




MeV/c 2 


Ttot ('72 (1950)") 


44111 


237 ±4 


587lf 




MeV 


r 77 £(7r 7r°) 




7fi+ 48 


3341™ 


(fixed) 


eV 


X 2 (ndf) 


323.2 (311) 


594.4 (313) 


1397.8 (315) 


2306.8 (315) 





the background amplitudes, and uncertainties from the unfolding procedure. 

For each study, a fit is made allowing all the parameters to vary. The differences of the 
fitted parameters from the nominal values are quoted as systematic errors. Again, 1000 sets 
of randomly generated input parameters are prepared for each study and fitted to search for 
the true minimum and for possible multiple solutions. Unique solutions are found repeatedly 
in all the cases. Once the solutions are found, several tens of repeated minimizations are 
needed to obtain fits that converge. With many parameters (24 — 26 here) to be fitted, the 
approach to the minimum is rather slow. 

Two fitting regions are tried: a higher one (1.74 GeV < W < 2.60 GeV) and a lower one 
(1.66 GeV < W < 2.40 GeV). The normalization error studies are divided into those from 
uncertainties of the overall normalization and those from distortion of the spectra in both 
| cos 8* | and W. For the overall normalization errors, fits are made with differential cross 
sections multiplied by (l±a e (w^ cosS *^), where cr e is the relative efficiency error. For distortion 
studies, ±4% (±3%) errors are assigned over the |cos#*| (W) range and differential cross 
sections are distorted by multiplying by 1 ± 0.1| cos6>*| =p 0.04 (1 ± 0.075IU =F 0.1575). 

For studies of the background (BG) amplitudes, each of the waves is changed to a first- or 
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a third-order polynomial except for the Go wave, where a first-order polynomial is introduced 
for both the real and imaginary parts of the amplitude. Parametrization uncertainties due 
to the phase convention where the Do and G2 background amplitudes are taken to be real 
are estimated by making the S and D 2 real instead and by introducing imaginary parts for 
the D and D 2 terms. 



TABLE VI: Systematic errors 





/ 4 (2050) 


72(1950)" 


Source 


Mass 


Itot 


r 77 s 7r o 7r o 


Mass 


Ttot 






(MeV/c 2 ) 


(MeV) 


(eV) 


(MeV/c 2 ) 


(MeV) 


(eV) 


W-r&nge 


+6 
—0 


+0 

—28 


+4.3 

—0 


+0 

— 14 


+0 

—87 


+0 

—27 


Normalization 


+0 
-2 


+0 
-1 


+1.0 
-0.8 


+2 
-0 


+5 
-0 


-1-1 1 
-3 


JDldD. | tUo (7 | 


+15 
-16 


+0 
-0.7 


+2.7 
-2.4 


+2 
-0 


+5 
-3 


+12 
-3 


Bias:W 




-1 


±1 


±0.2 


+3 
-2 


+2 
-0 


+ 1 
-0 


Unfolding 


+35 
-0 


+0 

-68 


+0 
-3.0 


+0 
-44 


+0 
-84 


+0 
-36 


BG: ReS 


+50 

-0 


+0 
-72 


+0 
-3.1 


+0 
-46 


+9 
-88 


+0 
-37 


BG: ImS 


+0 
-1 


+2 
-7 


+0 
-0.1 


±1 


+9 
-0 


+9 
-0 


BG: D 


+1 

-13 


+9 
-15 


+2.7 
-0 


+3 
-1 


+3 
-4 


+22 
-0 


BG: ReD 2 


+36 
-0 


+29 
-2 


+11.9 
-0.6 


+0 
-22 


+24 
-94 


+16 
-24 


BG: IrmD 2 


+0 
-12 


+0 
-12 


+3.1 
-0 


+0 
-4 


+0 
-13 


+11 
-0 


BG: Go 


+20 
-0 


+0 
-23 


+0 
-1.0 


+2 
-21 


+6 
-49 


+ 1 
-19 


BG: G 2 


+205 
-0 


+7 
-69 


+19.1 
-0 


±10 


+0 
-54 


+377 
-13 


BG: Real D & G 2 


±6 


+1 
-11 


±0.4 


+2 
-1 


+ 1 

-6 


+7 
-6 


Total 


+218 


+31 


+23.5 


+12 


+28 


+379 


-25 


-129 


-5.2 


-73 


-192 


-68 



The resulting systematic errors are summarized in Table IVIL Total systematic errors 
are calculated by adding the individual errors in quadrature. We obtain the mass, total 
width and r 77 i3(7r°7r°) of the / 4 (2050) to be 1884 ±\j tf* MeV/c 2 , 453 ± 20 MeV and 
7-7 ±U ±5 3 2 5 eV, where the first errors are statistical and the second systematic. The errors 
are dominated by systematics, and mostly come from uncertainties due to the unfolding 
procedure and background parametrization, and possible biases in the cos 8* distribution. 

From the measured branching fraction to titt (Table II Vj) . the two-photon width of the 
/ 4 (2050) is obtained to be 136 1| 2 1 91 5 eV. Given the large systematic error, we cannot 
conclude that the two-photon width of the / 4 (2050) is nonzero. However our data clearly 
require a G-wave component (see Fig. [6]), and the unacceptable fit without the / 4 (2050) 
(Table [V} strongly supports a finite two-photon coupling. In the past, TASSO and JADE 
have set 95% confidence upper limits for the / 4 (2050) to be T^B(KK) < 0.29 keV 
and r 77 £>(7T7r) < 1.1 keV 19 ]. respectively. These can be translated into upper limits for 



the two-photon widths of 43 and 6.5 keV, respectively. The power of such a large statistics 
(3 orders of magnitude more) of our experiment is evident. The nominal fit brings quite 
unexpected results: it requires a "flip" of the / 4 (2050) and "/2(1950)" positions with the 
mass of the former becoming 1885t 4 3 MeV or 153 MeV lower than the optimal mass of the 
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"/2(1950)". In addition, the fit requires both states to be much broader than before, 440-450 
MeV or about 2 times larger than their PDG values. Obviously, the interference of the D2 
and G2 amplitudes with each other and with the underlying continuum demands a more 
sophisticated description probably involving more than one resonance in each wave. Such a 
full amplitude analysis is beyond the scope of this work. On the other hand, results of all 
the fits provide unambiguous evidence for a nonzero two-photon coupling of the G wave. 



V. ANALYSIS OF THE HIGHER-ENERGY REGION 

In general, we expect that at high energies and large scattering angles, leading term QCD 
calculations give reasonable predictions for hard exclusive processes such as 77 — > MiM 2 , 
where Mu 2 ) are mesons. However, at what energies these terms begin dominating depends 
on the hadrons involved. In addition, even at the highest energies the differential cross 
sections depends on the shape of the Mi (2) wave functions. For charged meson pairs such 
as 7r + 7r~, K + K~, the differential cross sections is only slightly sensitive to the shape of 
the wave functions and the numerically largest term in the differential cross sections is 



proportional to sin" " 4 0* [9|, llOl . 111] ]. However for neutral meson pairs this term is absent; 
the cross section da/d\ cos 0* Mis much smaller and much more sensitive to the shape of the 
meson wave functions BUS, Ell- 
in contrast, the main idea of the handbag model 12j is that the terms that are asymp- 
totically power corrections, give the numerically largest contributions even at currently 
available energies. The universal prediction of the handbag model is that the ratios 
da(M°M°)/da(M + M~) are constant, i.e. the energy and angular dependences are the 
same for charged and neutral mesons. In particular, da(7T°TT°)/da(7r + 7r~) = 0.5 [12J while it 



varies from « 0.07 at cos0* = to « 0.04 at I cos 0*1 = 0.6 in Ref.llOj. 



A. Angular dependence 

We compare the angular dependence of the differential cross sections in the range 
I cos 0*| < 0.8 for W > 2.4 GeV with the function sin~ 4 0*. We also try a fit with an 
additional cos 2 0* term, to quantify a possible deviation from the sin -4 0* behavior. We 
choose this function because it gives relatively good fits over a wide range in W. Thus the 
fit function is parametrized as: 

da/d\cos9*\ = a(sin" 4 9* + bcos 2 9*). (9) 

We fit using a binned maximum likelihood method and 16 bins in the range | cos0*| < 0.8. 
We know that the effect of charmonia is large in the region 3.3 GeV < W < 3.6 GeV, but we 
cannot separate it in the angular dependence because we cannot assume here any functional 
shapes for the noncharmonium component. The results of the fit for b are shown in Fig. [HI 
as well as the fit to the angular distributions in the four selected W regions, where the 
differential cross sections, the vertical axis of this figure, are normalized to the total cross 
section a(\ cos0*| < 0.8) in each W region, i.e. the area under the curve is 1. The parameter 
b is close to zero above W > 3.1 GeV compared to b ~ 10, when the contribution of the 
b cos 2 9* term in the total cross section, <x(| cos0*| < 0.8), is comparable to the contribution 
of the sin -4 9* term. The b parameter becomes nearly constant and then systematically 
negative above the charmonium region. The change in the b parameter, which approaches 
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FIG. 8: (a) The fits of the angular dependence of the normalized differential cross sections (see 
text) at four selected W points. For the dashed curves the coefficient b (see the fit formula in the 
text) is fixed to 0. The solid curves show the fits with b floating, (b) The energy dependence of 
the parameter b giving the best fits. Here, the charmonium contributions are not subtracted, and 
the data in the XcO an d Xc2 charmonium regions are plotted with open squares. 



a constant value near zero, occurs at a W value close to that observed in the charged-pion 
case 



B. Yields of Xcj charmonia 

The structures seen in the yield distribution for 3.3 GeV < W < 3.6 GeV and | cos#*| < 
0.4 (Fig. [9]) are from charmonium production, 77 — > Xco, Xc2 —> 7r°7r°. Similar production of 
the two charmonium states is observed in the 7r + 7r~, K + K~ and K^Kg final states 0, H[. 

We fit the distribution to contributions from the XcCb Xc2 and a smooth continuum com- 
ponent using the following function: 

Y{W) = Iv^^+e^y^BW^^ (10) 
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in the W region between 2.8 and 4.0 GeV, where BW XcJ (W) is a Breit-Wigner function 
for the charmonium amplitude, which is proportional to 1/(W 2 — M Xcj — iM XcJ F XcJ ) and is 
normalized as / |BW XcJ (iy)| 2 (iiy = 1. The masses and widths, M and T, of the charmonium 
states are fixed to the PDG world averages 17| . The component aW~@ corresponds to the 
contribution from the continuum, with a fraction k that interferes with the Xco amplitude 
with a relative phase angle, <fi. It is impossible to determine the interference parameters for 
the Xc2 because of its much smaller intrinsic width compared to experimental resolution. We 
fit the Xc2 yield {N XC 2) with a formula where no interference term is included, and later we 
estimate the maximum effects from the interference term when determining the two-photon 
decay width of \ci- We use data only in the range | cos^*! < 0.4 where the charmonium 
contribution is dominant. Smearing effects due to a finite mass resolution are taken into 
account in the fit, using the same function as used for the unfolding. 

A binned maximum likelihood method is applied. We examined two cases with and 
without the interference. Reasonably good fits are obtained for both cases. The fit results 
are summarized in Table IVII1 In the table, L is the likelihood value and ndf is the number 
of degrees of freedom. The normalization iV xc o in Eq.(8) is proportional to the square of 
the resonance amplitude. The yields from the fits are translated into products of the two- 
photon decay width and the branching fraction, T^(xcj)^(Xcj vr°7r°), which are listed in 
Table IVIIII The systematic errors are taken from the changes in the central values of the 
fitted yields when the absolute energy scale is varied by ±10 MeV for the W measurement, 
the invariant-mass resolution is varied by ±10% for the corresponding Gaussian widths, and 
the fitting range is narrowed to the range 2.96 - 3.84 GeV, and when the efficiencies are 
varied by their uncertainties. The changes in the goodness of fit (— 21n£) for the first two 
variations are found to be small, less than 1.7, for the interference case. 

The Xco is observed with a statistical significance of 7.6a (7.3cr) when we take (do not 
take) interference into account. The statistical significance for the Xc2 is 2.6a when we take 
interference of the Xco m to account, but it is only 1.3cr when we do not take into account 
interference. This is because interference makes the line shape of the Xco highly asymmetric 
with a short tail and destructive interference on the high-energy side. The solid and dashed 
curves in Fig. [9] show the fits for the two cases (with and without Xco interference). 

The results for r 77 £>(x c j) in the 7r°7r° final state can be compared to the only direct 
measurement of this quantity in the 7r + n~ decay mode from Belle, 15.1 ± 2.1 ± 2.3 eV and 
0.76±0.14±0.11 eV for the Xco an d Xc2, respectively [2]. Although the effects of interference 
were neglected in the 7r + 7r~ measurements, the results are consistent with the ratio expected 
from isospin invariance, B(xcj — ► 7r ° 7r °)/'B(XcJ —> 7r + 7r~) = 1:2. Our results for T 77 i3(x c j) 
for the Xco(2) agree within errors with the indirect determination of these quantities using 
the corresponding world averages [l7| or recent measurements of £>(x c o(2) — ► 7r°vr°) [loj as 
well as B{xc.j 77) and r 77 (x c j) [2l[ by the CLEO collaboration. 
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TABLE VII: Results of the fits (see text) to obtain the charmonium contributions with and with- 
out interference effects. Errors are statistical only. Logarithmic likelihood (ln£) values are only 
meaningful when comparing two or more fits. 



Interference 


N v n k 

XcO 


(f> N xc2 


-2ln C/ndf 


Without 
With 


100 ± 16 

103121 0.82t°i* 


(1.1±0.3)tt 34 ± 13 


52.4/56 
44.2/54 



TABLE VIII: Products of the two-photon decay width and the branching fraction for the two 
charmonia. Here, T J1 B(xcj) means r 77 (x c j)£>(Xcj — > 7r°7r°). The first, second and third errors 
(when given) are statistical, systematic and from the maximal uncertainties of the relative phase 
in Xc2 production. 



Interference 


rV?( Xc0 ) (eV) r 77 £( Xc2 ) (eV) 


Without 
With 


9.7 ±1.5 ±1.2 0.18±g;i| ± 0.08 
9.9±|;g ± 1.6 0.48 ±0.18 ±0.07 ±0.14 



C. Subtraction of the charmonium contributions 

We subtract the charmonium contributions from nearby bins of the charmonium (Xcj) 
region, 3.3 — 3.6 GeV, in order to obtain a pure differential cross section from the continuum 
component. We use the fit result with interference obtained in the previous subsection. 

The estimated charmonium yield that includes the contribution from the interference term 




FIG. 9: The W distribution of the candidate events with | cos 8*\ < 0.4 near the charmonium region. 
The solid and dashed curves show the fits described in the text with and without interference with 
the XcO- 
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TABLE IX: The value n in c tot oc W n in various reactions fitted in the W and | cos#*| ranges 
indicated. 



Process 


n 




W range (GeV) 




cos 9*\ ranj 


^e Reference 


vr 7r° 


6.9 ±0.6 ±0.7 


3.1 


-4.1 (exclude 3.3 


- 3.6) 


< 0.6 


This exp't 


7T + 7T~ 


7.9 ± 0.4 ± 1.5 




3.0 - 4.1 




< 0.6 


m 


K+K- 


7.3 ± 0.3 ± 1.5 




3.0 - 4.1 




< 0.6 


[2] 




10.5 ±0.6 ±0.5 


2.4 - 


4.0 (exclude 3.3 - 


3.6) 


< 0.6 


\£ 


vr 7r° 


8.0 ±0.5 ±0.4 


3.1 - 


4.1 (exclude 3.3 - 


3.6) 


< 0.8 


This exp't 



is converted to a differential cross section contribution in each angular bin of | cos6**| < 0.8 
by assuming a flat distribution for the Xco component and a distribution ~ sin 4 9* for the 
Xc2 component This assumption is only a model. In fact, we do not know the angular 
distribution of the interference term; the charmonium amplitudes can interfere with the 
continuum components with different J's of unknown sizes. 

For the W = 3.25 GeV bin, the fit result indicates that there is a non-negligible effect 
from the Xco when we assume interference, and thus we make a correction for charmonium 
subtraction. The contribution of the charmonium components in the original differential 
cross sections is 18% at | cos#*| < 0.6. For W = 3.3 — 3.6 GeV, we apply a subtraction for 
the angular bins 0.4 < | cos#*| < 0.8 after extrapolating the charmonium yield determined 
in the range | cos#*| < 0.4. 

The differential cross section thus obtained for the continuum is integrated over the range 
| cos9*\ < 0.6. We convert a (OA < | cos9*\ < 0.8) to a(\ cos9*\ < 0.4) for W = 3.2-3.6 GeV, 
by assuming that the angular dependence of the differential cross section is ~ sin -4 9*. The 
results are plotted in Fig. [TOTa). where the cross section for 77 — > 7T + n~ from Ref. 0] is also 
shown. 

D. W dependence and ratio of cross sections of 7r°7r to tt + tt~ 

We fit the differential cross sections integrated over the polar angle, <r(| cos#*| < 0.6), to a 
power law in the cm. energy, W~ n , for the energy region 3.1 GeV < W < 4.1 GeV, in which 
the angular dependence of the differential cross section does not show any large changes. 
In the fit, we do not use the data in the charmonium region (W = 3.3 — 3.6 GeV), where 
we cannot determine the cross section of the continuum component in a model-independent 
manner. 

The result of the power-law fit is summarized in Table IIXI and compared to that for 
other processes. The systematic error is dominated by the uncertainty of the charmonium 
contribution in the range 3.1 GeV < W < 3.3 GeV. This n value is compatible with the 
results for the 7r + 7r~ and K + K~ processes [2j], but significantly smaller than that in the 
K° S K° S case fl. 



The fit for 3.1 GeV < W < 4.1 GeV is shown in Fig. [POT a). which also shows the 
cross section of 7r°7r° from the Belle measurement [2]. In Fig. [TOT b) we show the ratio of 
the cross sections of 7r°7r° to ir + 7r~. This ratio is rapidly falling at low energies, but its 
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2.2 2.7 3.2 3.7 4.2 



W (GeV) 

FIG. 10: (a) The cross sections for the 77 — > 7r°7r° (solid circles) and 77 — > tt + it~ (triangles, 0]) 
for I cos 0*| < 0.6. The curve is a fit to the cross section for 77 — ► 7r°7r° with a ~ W~ n functional 
shape, (b) Ratio of the cross section for the tt°tt° process to the -k + -k~ process. The error bars 
are statistical only. The shorter horizontal line is the average for 3.1 GeV < W < 4.1 GeV. 
The horizontal line (0.5) is an expectation from isospin invariance for a pure 1 = component. 
In (a) and (b), the estimated charmonium contributions are subtracted in both it + it~ and 7r°7r° 
measurements. The results in the W region 3.3 - 3.6 GeV (plotted with gray circles) are not used 
for the fits. 

behavior changes above 3.1 GeV, where the two processes have similar W~ n dependence, 
which results in the almost constant ratio. The average of the ratio in this energy region is 
0.32 ± 0.03 ± 0.05, where the data in the 3.3 - 3.6 GeV region is not used when calculating 
this average. This ratio is significantly larger than the prediction of the leading-order QCD 
calculations (9I Hoi fill and is somewhat smaller than the value of 0.5, which is suggested by 
isospin invariance [121 ]. 
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VI. SUMMARY AND CONCLUSION 



We have measured the process 77 — > 7r°7r° using a high-statistics data sample from e + e~ 
collisions corresponding to an integrated luminosity of 223 fb -1 collected with the Belle 
detector at the KEKB accelerator. We derive results for the differential cross sections in the 
center-of-mass energy and polar angle ranges, 0.6 GeV < W < 4.1 GeV and | cos^*! < 0.8. 

Differential cross sections are fitted in the energy region 1.7 GeV < W < 2.5 GeV in a 
model where the partial waves consist of resonances and smooth backgrounds. Various fits 
are performed that provide unambiguous evidence for a nonzero two-photon coupling of the 
G wave. Helicity-2 production (G2) is preferred compared to the helicity-0 (Go) one. 

We observe production of the charmonium state Xco and obtain the product of its two- 
photon decay width and the branching fraction to 7r°7r°. The angular distribution of the 
differential cross section is largely energy dependent, and approaches ~ sin _4 #* above W = 
3.1 GeV. This observation and the energy dependence of the cross section above this energy, 
which is well fitted by W~ n , n = 6.9 ± 0.6 ± 0.7, are compatible with those measured in the 
7T + 7r~ channel. We obtain the cross section ratio, o"(7r 7r°)/cr(7r + 7r~), to be 0.32 ± 0.03 ± 0.05 
on average in the 3.1-4.1 GeV region. This ratio is significantly larger than the prediction 
of the leading-order QCD calculation. 
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